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Cesium adsorption structures on Ag(lll) were characterized in a low-temperature scanning tun- 
neling microscopy experiment. At low coverages, atomic resolution of individual Cs atoms is occa- 
sionally suppressed in regions of an otherwise hexagonally ordered adsorbate film on terraces. Close 
to step edges Cs atoms appear as elongated protrusions along the step edge direction. At higher 
coverages, Cs superstructures with atomically resolved hexagonal lattices are observed. Kinetic 
Monte Carlo simulations model the observed adsorbate structures on a qualitative level. 
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I. INTRODUCTION 

Alkali metal adsorption on metal surfaces has been 
studied intensely since many years (for a review the 
reader is referred to Refs. 1,2). With their single s elec- 
tron at the outermost atomic shell these metals are con- 
sidered as simple and much research is devoted to the 
analysis of the interplay between their mutual interaction 
and the interaction with the hosting substrate.^ Previ- 
ous structural studies of Ag(lll)-Cs focussed on adsor- 
bate superstructures at higher coverages. For instance, 
a p{2 X 2) structure along with evidence of some disor- 
dered and possibly incommensurate phases was found. 
Bond length changes between the Cs adsorption layer 
(adlayer) and the Ag(lll) substrate surface were studied 
by surface-extended X-ray absorption fine structure also 
for higher coverages.^ Leather man and Diehl provided a 
structural analysis of Ag(lll)-Cs for various coverages 
and temperatures using low-energy electron diffraction.^ 
In particular, for very low coverages they found ring- 
like diffraction patterns which were assigned to disor- 
dered or fluid over layer phases. This fluid phase appears 
to be a common arrangement of alkali metal films on 
metal surfaces at low coverage. ''^^^^^^^ It is usually ar- 
gued that at low coverages the dipole-dipole repulsion 
between the adsorbed alkali metal atoms dominates the 
alkali-substrate interaction resulting in a structure with 
no long-range order, yet with a typical distance between 
the adsorbed atoms (adatoms). Besides these structural 
properties, alkali metal-covered surfaces exhibit intrigu- 
ing electronic properties resulting from quantum size ef- 
fects: Quantum well states confined to ultrathin films of 

alkali metals are investigated experimentally as well as 
theoretically^a^a^^a^a^Oia^aa^^^^^ 

While the dipole-dipole interaction presents a direct 
coupling between Cs adatoms, the substrate may influ- 
ence the Cs-Cs interaction in an indirect way. Indirect in- 
teraction between two atoms was first investigated theo- 
retically by Koutecky^ and then by Grimley,^^ Newns,^^ 



and Einstein and Schrieffer.^^ In an early field ion mi- 
croscopy experiment evidence for indirect interaction be- 
tween Re atoms adsorbed on W(llO) was reported by 
Tsong.^^ In particular, it was found that the interac- 
tion energy exhibits an oscillatory behavior as a function 
of the separation distance between the adatoms. Lau 
and Kohr^ then predicted a long-range and oscillatory 
interaction between atoms mediated by Friedel oscilla- 
tions of a two-dimensional electron gas. The first scan- 
ning tunneling microscopy (STM) experiment evidenc- 
ing this type of interaction was reported by Brune et 
al.^^ for carbon atoms adsorbed on Al(lll). Similar 
experiments were then performed for benzene molecules 
on Cu(lll),^^ sulfur atoms on Cu(lll),^^ copper atoms 
on Cu(lll),33 ^^-^^Yi atoms on Cu(lll) and Ag(lll)^ 
and for Ce atoms adsorbed on Ag(lll).^^ These exper- 
iments reveal that the adatom-adatom distance is pre- 
dominantly influenced by the surface state-mediated in- 
teraction. In particular, the mutual distance is reported 
to be Af/2 where Ap denotes the Fermi wavelength of 
the involved surface state. This interaction may lead 
to superlattices on surfaces^ or to confinement- induced 
adatom self-organization in quantum corrals ^^S*^ 

Here, we report a low-temperature STM experiment 
on Cs adsorbed on Ag(lll) at various coverages. At 
low coverages, substrate terraces are covered by Cs 
adatoms which locally exhibit hexagonally ordered do- 
mains. Long-range hexagonal order of the adlayer, how- 
ever, does not exist owing to regions where individual 
Cs adatoms are not resolved in STM images. Similarly, 
at and close to step edges individual Cs atoms are not 
resolved and appear as almost continuous rows running 
parallel to the step edges. At higher coverages the char- 
acteristic smearing of adsorbate structures on terraces 
and at step edges disappears and individually resolved 
adatoms form superlattices with long-range hexagonal 
order. Kinetic Monte Carlo simulations model the exper- 
imental situation on a qualitative level. We found that 
both the surface state-mediated interaction between Cs 
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adatoms and the dipole-dipole interaction are equally im- 
portant, in particular for low coverages, and have to be 
included in the model. The simulations reveal that loss of 
atomic resolution is likely due to the low adsorbate cov- 
erage matching an incommensurate superstructure. At 
higher coverages, as a result of a reduced adatom-adatom 
distance, the mutual dipole repulsion is strong enough to 
stabilize a superlattice with long-range hexagonal sym- 
metry and to suppress adatom diffusion. 



II. EXPERIMENT 



with Aq denoting the scattering amplitude, £'0 the sur- 
face state binding energy, Sq the scattering phase shift, 
and the Fermi wave vector. 

Our idea is that regions of the adsorbate lattice that 
exhibit loss of atomic resolution are due to an enhanced 
mobility of the adatom. It is thus desirable to access 
the probability of finding a particle at a given lateral co- 
ordinate. Therefore, for comparison with experimental 
data we use the time-averaged two-dimensional distribu- 
tion function, ^r(r), which represents the probability of 
finding a particle at a position r during a time interval 



Experiments were performed using a custom-built 
scanning tunneling microscope operated at 7K and in 
ultrahigh vacuum with a base pressure of 10~^ Pa. The 
Ag(lll) surface as well as chemically etched tungsten 
tips were cleaned by argon ion bombardment and anneal- 
ing. Cesium was deposited at room temperature from 
commercial dispensers^^ at a rate of ^ 0.05MLmin~^ 
as monitored by a quartz microbalance and judged from 
the deposition time and corresponding STM images. We 
define a monolayer (ML) as one Cs atom per Ag atom. 
All STM images presented in this work were obtained in 
the constant-current mode with voltages applied to the 
sample. 



III. SIMULATIONS 

We model aspects of our results by using the kinetic 
Monte Carlo method^^^^^ which got a firm theoretical 
background by Fichthorn and Weinberg. Recently, the 
kinetic Monte Carlo method has been used to study 
the self-organization of adatoms due to surface state- 
mediated interactions i^i^ The substrate lattice is mod- 
eled by finite hop rates of adatoms between adjacent lat- 
tice sites ri and rj. The hop rate Uij is described by an 
Arrhenius law according to 



uo ex.p{-Eij/kBT) , 



(1) 



where T is the substrate temperature, ke Boltzmann's 
constant and uq the attempt frequency. For the hop- 
ping barrier between lattice sites i and j we use Eij = 
{Ej—Ei) where is the diffusion barrier height 
for an isolated adatom and the £j(i) describe the inter- 
action of atom j (i) with all the other atoms. The total 
interaction is the sum of the dipole-dipole coupling 
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where p is the dipole moment of an individual Cs adatom, 
60 the dielectric constant, and the surface state-mediated 
interaction^ 



El = A^E^ 
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where N is the number of particles used in the simulation. 



IV. RESULTS AND DISCUSSION 

A. Low coverage: < 0.1 ML 

Figure [1] illustrates the main experimental findings at 
low coverages. The STM image in Fig.llfa) shows two 
adjacent terraces of a Ag(lll) surface covered with 0.03 
- 0.04 ML Cs. On terraces we observe local hexago- 
nal order of Cs atoms. Long-range order is suppressed 
in regions where Cs adatoms are no longer resolved as 
point-like features [Fig.fljb)]. Rather, these regions are 
characterized by broadened protrusions extending along 
the symmetry directions of the adsorbate lattice. From 
atomically resolved STM images of the adsorbate lattice 
a mutual adatom distance of (1.5 ± 0.2) nm was deter- 
mined. At and close to step edges [Fig.[TJc)] adatoms 
tend to form rows parallel to the step direction sepa- 
rated by (1.5 ± 0.2) nm. Individual Cs adatoms start to 
be resolved again at distances exceeding 5 nm from the 
step edge. 

Below we suggest that the peculiar coexistence of re- 
gions with hexagonal order together with regions in which 
adatoms are not resolved is likely due to an adsorbate 
coverage matching an incommensurate phase. For our 
kinetic Monte Carlo simulations several parameters have 
to be estimated, which are discussed in the following. 
The results of the kinetic Monte Carlo simulations de- 
pend on the transition rates, Eq. ([T]), which are directly 
coupled to the parameters of the dipole and the surface 
state-mediated interactions. These parameters have to 
be estimated in advance. 

In a first step we estimated the dipole moment of an in- 
dividual Cs adatom. Adsorption leads to a considerable 
charge transfer from an alkali atom to the substrate. As 
a consequence the adsorbed atom becomes partly ionic 
and develops a dipole moment which leads to a decrease 
of the work function. From the initial linear decrease 
of the work function the dipole moment of the adsorbed 
species may be evaluated^ The apparent height of the 
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FIG. 1: (Color online) (a) STM image of two adjacent ter- 
races of Ag(lll) covered with 0.03 - 0.04 ML Cs deposited 
at room temperature and imaged at 7K. A monatomically 
high step separates the terraces (voltage V = 200 mV, cur- 
rent / = 0.2 nA, size 97 nm x 97 nm). (b) Close-up view of 
Cs adatoms on a terrace (35 nm x 35 nm). (c) Close-up view 
of the adsorbate arrangement in the vicinity of a step edge 
(20 nm X 20 nm). 

tunneling barrier in STM being related to the work func- 
tion, we acquired current-distance curves on clean and 
Cs-covered Ag(lll) from which we inferred apparent bar- 
rier heights of 5 eV and ^ 4 eV, respectively. The 1 eV 
difference of the apparent barrier heights is in good agree- 
ment with the work function modifications observed on 
other surfaces covered with alkali metals. ^^^^^^^ Using a 
work function change of A$ ^ 1 eV we extract a dipole 
moment p according tc^^ 

eoA<l> = epn, (5) 

where — e is the electron charge and n the surface density 
of alkali metal atoms. As a result we obtain p (0.09 ± 
0.03) enm= (4.3 ± 0.7) D. 

To compare the strength of the dipole-dipole inter- 
action [Eq. ([2])] with the interaction mediated by the 
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FIG. 2: Dipole-dipole interaction energy (gray) and sur- 
face state-mediated interaction energy (black) between Cs 
adatoms on Ag(lll) calculated according to Eqs. ([2|) and 
(|3]) with parameters p — 0.09 enm, Aq — 0.3, — 7r/2, 
Eq = — 0.06eV, and /cf = 0.813 nm~^. At the experimentally 
observed lattice constant of the adlayer, a ^ 1.5 nm (dashed 
line), the surface state- mediated interaction and the dipole- 
dipole interaction energy have similar values. 



surface state [Eq. (|3|)] we also have to determine the 
parameters Aq, Eq, and ^o- From previously pub- 
lished results^^'^i^i^ we infer an amplitude of Aq ~ 
0.3. The scattering phase shift varies between 0.37r and 
Q^^ ^32^33^34^35 within this range our simulations do 
not depend strongly on the specific choice of ^o- We 
therefore decided to take 5o = 7r/2. For clean Ag(lll) 
the surface state binding energy is Eq ~ — 0.06eV below 
the Fermi energy. Charge transfer from the Cs layer to 
the substrate, most likely induces a shift of the surface 
state binding energy to higher values. To estimate this 
shift, we resort to a similar adsorption system, namely 
Cu(lll)-Cs, for which a change of the Cu(lll) surface 
state binding energy has been investigated by photo- 
electron spectroscopy.^^ By extrapolating the photoemis- 
sion results to low coverages, we find an energy shift of 
~ 0.01 eV to higher binding energies for a coverage of 
0.03 ML. However, taking Eq ^ -0.07 eV rather than 
£^0 ^ — 0.06eV does not alter the results of our simula- 
tions. For the parameters estimated above we calculated 
the interaction energies for the dipole-dipole and the sur- 
face state-mediated coupling according to Eqs. ([2|) and 
([3|). A comparison is shown in Fig.O At the adsorbate 
lattice spacing of a ^ 1.5 nm observed in the experiment 
(dashed line in Fig. [2]) both interaction energies are sim- 
ilar. Therefore, in our simulations of the low-coverage 
adatom arrangements we take both interactions into ac- 
count. 

To estimate the diffusion barrier height, £^d, we ob- 
served single Cs atoms on the cold surface. Upon po- 
sitioning the tip above the center of a Cs adatom the 
feedback loop of the instrument was opened. Next, we 
measured the time interval from the moment of freez- 
ing the feedback loop until the current dropped, which 
signals that the adatom has moved away from its orig- 
inal position. This procedure resulted in a distribution 
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FIG. 3: (Color online) Density plots of the time- averaged 
distribution function Qrir) for T = 7K, = 17meV, and 
T = 12 s calculated for particle coverages of (a) 0.037 ML and 

(b) 0.032 ML. The image size corresponds to 30 nm x 30 nm. 

(c) Time evolution of the fraction of defects, i. e., of parti- 
cles with 5 and 7 nearest neighbors for an incommensurate 
(circles) and a commensurate (squares) phase. 



of time intervals between 1 and 15 s. Together with an 



assumed attempt frequency of uq 
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we estimate a diffusion barrier height between 16 and 
18meV at 7K. We find that the choice of = 10^^ s"^ 
and ^ 17meV matches the time scale of the experi- 
ment quite well. However, a word of caution is necessary 
in this context. Simulated time scales vary with the fac- 
tor i/Q exp(— £^D/kB^)- Consequently, by modifying vq^ 
Ed, or both within reasonable limits, experimentally ob- 
served time scales may be reproduced. If experimentally 
observed time scales differ from numerically simulated 
ones by a factor of more than 10^ then it is no longer 
reasonable to compensate this discrepancy by an adjust- 
ment of uo and £^d-^^ In this case new important physi- 
cal processes have to be included in the menu of hopping 
events of the kinetic Monte Carlo algorithm4i It turns 
out that the inclusion of the adatom mutual interactions 
according to Eqs. ([2]) and (|3]) already leads to correlated 
transition processes. 

With these parameters at hand we then addressed the 
experimentally observed coexistence of ordered and less 
ordered regions. To monitor the time evolution of an 
initial particle configuration we calculated first the ther- 



mal equilibrium of a particle lattice at 6 K using canon- 
ical Monte Carlo simulations. In the subsequent kinetic 
Monte Carlo calculations we set the temperature to 7K 
as measured in the experiment and analysed the time 
evolution of the particle lattice. Figures [31(a) and E^b) 
show the results for particle coverages matching a com- 
mensurate [0.037 ML, {V27 x ^27)^^30°] and an incom- 
mensurate (0.032 ML) superlattice, respectively. The fig- 
ures show plots of ^r(r) with r = 12 s. Surprisingly, such 
a slight variation in the coverage leads to considerably 
different results. Whilst in the commensurate adsorp- 
tion phase [Fig.[3Ka)] gr{'^) is characterized by localized 
and sharp maxima at sites of a hexagonal lattice with 
long-range order, in the incommensurate phase [Fig.lSJb)] 
^r(i*) exhibits regions with maxima which are broader 
than in the commensurate phase. In particular, at some 
positions adjacent maxima overlap and appear as elon- 
gated and continuous lines following close-packed direc- 
tions of the hexagonal particle lattice. As a first result 
we may summarize that particle arrangements that are 
incommensurate exhibit enhanced particle mobility. 

We performed a Voronoi analysis^*^ to investigate 
the relation between ^r(r) and particle mobility in more 
detail. As a result we obtain the fraction of particles 
with a number of nearest neighbors deviating from six. 
These particles are referred to as defects in the following. 
This analysis thus provides access to local distortions of 
the hexagonal symmetry. We define the fraction of these 
particles according to Qr = N^/irN) with denoting 
the number of defects recorded during the time interval 
r. Figure [3Kc) shows the time evolution of qr for the 
commensurate (squares) and the incommensurate (cir- 
cles) phase shown in Figs. [31(a) and [3Kb). The fraction 
increase for < r < 6s is related to thermal relaxation 
of the initial particle arrangement at low temperature. 
The relaxation period is finished at r 12 s and qr then 
stays essentially constant for the commensurate phase. 
For the incommensurate phase, however, an increase of 
qr is still observed until r = 25 s. Due to costly com- 
puting time we have not performed kinetic Monte Carlo 
simulations for larger r but we expect qr to become con- 
stant also for the incommensurate phase reflecting then 
thermal equilibrium of the superstructure. From these 
simulations we may conclude that once the fraction qr 
does not exceed 5 - 10 % the superlattice remains stable 
[Fig. [31(a)] while larger values of qr indicate the onset of 
lattice melting [Fig. [31(b)] or thermal equilibrium which 
has not been reached. 

To model the adsorbate arrangement at step edges, the 
steps are taken into account as scatterers of the surface 
state. We considered a step edge as an infinitely long 
and linear chain of Ag atoms each of which gives rise 
to a surface state-mediated interaction with Cs adatoms 
according to Eq. (j3|). Integrating the contribution of each 
scatterer leads to the total interaction^^' 
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FIG. 4: (Color online) Density plot of 5'T(r) in the vicinity of 
a (111) (upper part) and a (100) (lower part) step edge. The 
simulation area corresponds to an image size of 30 nm x 30 nm. 
The density plot was generated for T = 7K, r = 12s, and 
Ey^ = 17meV. 



where r±_ denotes the distance to the step and ao = 
0.287 nm the Ag(lll) lattice constant. The parameters 
Eq^ /cf, and 5q are the same that we used for model- 
ing the adlayer structure on terraces, while the choice of 
^0 = 4meV is in good agreement with the results ob- 
tained for Cu atoms adsorbed on Cu(lll)^ Figure [H 
shows the density plot of ^r(r) for r = 12 s. Two kinds 
of step edges were included in the simulations, namely 
a step edge of the (111) type (upper part of Fig.[4|) and 
a (100) step edge (lower part of Fig.[4|). Deviating from 
our simulations of adsorbate structures on terraces we 
found that the surface state-mediated interaction accord- 
ing to Eq. (|6|) introduced an additional stabilization of the 
super lattice. The distribution function, ^r(r), exhibits 
overlapping maxima along the step edge directions which 
is in agreement with Cs adatoms showing elongated rows 
close to step edges in STM images [Fig.[Hc)]. Addition- 
ally, maxima of ^r(r) become more localized again with 
larger distances from the step edge and reflect then the 
characteristic particle arrangement on terraces. 

According to the Monte Carlo simulations we interpret 
the experimental observations as a consequence of an en- 
hanced Cs adatom mobility at the low-coverage regime 
discussed in this section. We can identify the origin of the 
observed mobility as a property of coverages that match 
incommensurate adsorbate phases. 



B. Higher coverages: > 0.1 ML 

A way to stabilize the adsorption lattice, i. e., to ob- 
tain a hexagonal superstructure with long-range order, is 
the increase of the coverage. Many structural analyses of 
alkali metal films adsorbed on metal surfaces have pre- 
viously reported the transition from a disordered phase 
at very low coverage to well-ordered adsorbate structures 
at higher coverages (Ref.|l| and references therein). An 
increased coverage leads to a smaller adatom-adatom dis- 
tance and will therefore increase the dipole-dipole repul- 
sion between the adatoms. As a consequence, the proba- 
bility for a Cs atom to hop from one adsorption site to an 
adjacent one is lowered. Our experiment and calculations 
corroborate this picture. 

Figure [5^a) shows an STM image of Ag(lll) covered 
with 0.11 ML of Cs. We observe a hexagonal adlayer 
with long-range order and a mutual adatom distance of 
(0.74 ± 0.04) nm. This adsorbate superstructure corre- 
sponds to a (3 X 3) commensurate phase and has been 
reported before by low-energy electron diffraction.^ The 
inset of Fig. [51(a) shows the atomically resolved Ag(lll) 
lattice to indicate that the adlayer and substrate lattice 
exhibit the same orientation [see dashed line in Fig.[5](a)]. 
Regions indicating an enhanced mobility of adatoms are 
no longer present at this coverage. Given that the com- 
mensurate adsorbate phase is observed experimentally 
at a specific coverage, we performed simulations for an 
adatom lattice matching the substrate lattice. The super- 
structure similar to the one observed in the experiments 
was obtained as the ground state from canonical Monte 
Carlo calculations. The distribution function ^r(r) for 
this coverage was then calculated for a variety of time 
intervals. Even for extended intervals, e. r = 35 s, 
no enhancement of the diffusion processes was observed 
indicating that the increased dipole-dipole coupling sta- 
bilizes the hexagonal adsorption lattice. 

The same behavior was observed at higher coverages. 
FigurelSJb) shows an STM image of Ag(lll) covered with 
0.15 ML Cs. Again, no indication of enhanced diffusion 
of Cs adatoms was found. The Cs adatoms are arranged 
in an incommensurate hexagonal superlattice with a mu- 
tual distance of ~ 0.66 nm and rotated with respect to 
the Ag(lll) lattice by ~ 19°. Moreover, at this coverage 
domains of hexagonally ordered Cs atoms were observed. 
Cesium atomic rows of adjacent domains are offset by half 
an adlayer lattice constant [see inset of Fig.[5^b)]. At the 
domain boundaries rows of irregulary shaped structures 
occur whose apparent height depends on the tunneling 
voltage. The formation of domains may be understood 
in terms of stress release at the domain boundaries. Sim- 
ilar observations of translational adsorbate domains have 
been reported for oxygen adsorption on W(llO)^^ and on 
Rh(lll).^^ 

To model these observations within a kinetic Monte 
Carlo approach we simulated particle arrangements for 
coverages ranging between 0.11 ML and 0.2 ML (Fig. [6]). 
After an interval of r = 12 s is elapsed the coordinates 
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FIG. 5: (Color online) (a) Quasi-three-dimensional represen- 
tation of constant-current STM image of Ag(lll) covered 
with 0.11 ML of Cs (y = 0.25 V, / = 0.1nA, 14nmx 14 nm). 
Cesium adatoms which appear as almost circular protrusions 
exhibit a mutual distance of 0.74 nm. Inset: Atomically re- 
solved Ag(lll) lattice. The dashed line indicates that adlayer 
and substrate lattice have the same orientation, (b) STM im- 
age of Cs-covered Ag(lll) at 0.15 ML {V = 1.2 V, / = 0.2 nA, 
80.4 nm x 80.4 nm). Lines on terraces which contain irregu- 
larly shaped structures appearing as depressions at the ap- 
pUed tunneling voltage are boundaries between translational 
domains of the adsorbate lattice. Inset: Atomically resolved 
Cs domains. While the orientation of the adjacent domains 
is identical, the lattices are translated by half a superlattice 
constant (see dashed line). 



FIG. 6: Arrangement of 340 particles according to Kinetic 
Monte Carlo simulations using a time averaging interval of 
r 12s and a substrate temperature of T = 7K. The cover- 
ages are indicated at the top right of each plot. The horizontal 
(inclined) arrow indicates the crystallographic direction of the 
substrate (particle lattice). The rotation angle is defined as 
the smaller angle enclosed by the two arrows. 



of 340 particles with a coverage of 0.11 ML are given by 
the upper left plot in Fig.[6l This coverage corresponds 
to the (3 X 3) commensurate adsorption phase and we 
see that the superlattice aligns with the crystallographic 
direction of the substrate lattice (depicted as the hori- 
zontal arrow). With increasing coverage we observe that 
this alignment weakens. At 0.132 ML some domains are 
still oriented along the substrate crystallographic direc- 
tion while others enclose an angle of ~ 20° indicating that 
at 0.132 ML commensurate and incommensurate phases 
coexist. This coexistence is in agreement with previ- 
ous reports based on low-energy electron diffraction.^ At 
0.136 ML all particle domains exhibit a rotation angle of 
^ 20°. This angle is maintained until at 0.2 ML the par- 
ticle lattice is aligned with the substrate lattice again. In 
particular, the incommensurate phase observed experi- 
mentally at 0.15 ML with a rotation angle of ^ 19° is 
well reproduced by our simulations. We summarize cal- 
culated and measured rotation angles as a function of 
the coverage in Fig.H We notice that our simulations re- 
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FIG. 7: (Color online) Rotation angle of the superlattice with 
respect to the substrate lattice. Experimental data are de- 
picted as circles and calculated data as triangles. The dashed 
Une indicates a coverage at which two phases with different 
rotation angles coexist. 



produce the experimentally observed rotation angles for 
commensurate adsorbate phases. The continuous rota- 
tion of incommensurate phases reported for Ag(lll)-Cs^ 
and Cu(lll)-Cs^ however, is not concerned here. 

Besides the rotation angle of commensurate super- 
structures the simulations indicate the origin for do- 
main formation as observed experimentally at 0.15 ML 
[Fig.[5](b)]. Increasing the particle density in the simu- 
lations to 0.14 ML and thus being close to the (a/7 x 
\^)R19.1° commensurate phase leads to a particle ar- 
rangement where the domains have almost disappeared 
in favor of a nearly homogeneous particle film with 
long-range hexagonal order. At a particle coverage of 
0.155 ML the boundaries of adjacent domains overlap 
leading to an increase of the particle interaction at the 
boundaries. A further increase of the particle density to 



0.2 ML induces a decrease of the domain size and leads to 
a similar particle arrangement as observed for 0.136 ML. 
As a result we obtain that particle coverages giving rise 
to incommensurate phases exhibit the propensity to form 
domains, while commensurate phases lead to homoge- 
neous films. 



V. SUMMARY 

Cesium films adsorbed on Ag(lll) with sub monolayer 
coverage have been investigated by scanning tunneling 
microscopy. The resulting adsorbate arrangements are 
modeled by kinetic Monte Carlo simulations. At very low 
coverages, long-range hexagonal order of the superstruc- 
ture does not exist while at higher coverages a stable and 
hexagonally ordered adsorbate lattice is observed. Mod- 
eling of the adsorbate arrangement at very low coverages 
demonstrates the importance of both dipole-dipole repul- 
sion and surface state-mediated interaction between Cs 
adatoms. Long-range hexagonal order is inhibited if the 
coverage matches an incommensurate adsorption phase 
or a disorder-to-order transition. At higher coverages, 
the dipole-dipole repulsion between the Cs adatoms be- 
comes large enough to stabilize hexagonal superlattices. 
While commensurate adsorption phases are character- 
ized by a homogeneous adsorbate films, incommensurate 
phases exhibit domain-like patterns. The experimentally 
observed rotation of commensurate adsorbate lattices is 
modeled by our calculations. 
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